Density-functional theory for fermions in the unitary regime 
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In the unitary regime, fermions interact strongly via two-body potentials that exhibit a zero range 
and a (negative) infinite scattering length. The energy density is proportional to the free Fermi gas 
with a proportionality constant f. We use a simple density functional parametrized by an effective 
mass and the universal constant £, and employ Kohn-Sham density-functional theory to obtain the 
parameters from fit to one exactly solvable two-body problem. This yields £ = 0.42 and a rather 
large effective mass. Our approach is checked by similar Kohn-Sham calculations for the exactly 
solvable Calogero model. 
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Ultracold atomic Fermi gases have received consid- 
erable experimental and theoretical interest since the 
achievement of Fermi degeneracy by DeMarco and Jin 
. One of the most interesting features is the ability to 
tune the inter-particle interaction itself via a Feshbach 
resonance, and to study the system as it evolves from 
the BCS regime with weakly attractive interaction to- 
ward the point where the interaction is strong enough 
to form di-atomic molecules, and the system may un- 
dergo Bose- Einstein condensation (BEC). Very recently, 
the BEC-BCS crossover has been subject of numerous 
experimentaljljl [1 MM 0, II H El and theoretical 
studies ElllMfQllir 

The dividing point of this BCS-BEC crossover defines 
the unitar y reg ime and is of particular interest |l& ITtI 
liH liH l20l l2l| . Here, the two-body system exhibits a 
bound state at zero energy, and the two-body scattering 
length diverges. For dilute systems, the inter-particle 
distance is much larger than the range of the interaction 
and much smaller than the scattering length. Thus, the 
inter-particle distance is the only relevant length scale, 
and the energy must be proportional to that of a free 
Fermi gas, 



E{N) =£E TF {N) 



(1) 



the dimensionless proportionality constant being denoted 
as £. Here, E(N) is the energy of the fully paired N- 
fermion system while Etf(N) is the Thomas-Fermi en- 
ergy of N noninteracting spin-1/2 fermions. The con- 
stant £ is universal, as it describes the physics of any 
dilute Fermi gas in the unitary regime. Approximately, 
it also describes systems close to the unitary regime, and 
can also be applied to dilute neutron gases as the two- 
neutron system also exhibits a scattering length that is 
much larger than the range of the nucleon-nucleon in- 
teraction. The exact determination of the universal con- 
stant £ is thus an important task. 

Recently, this constant was reliably determined 
through Monte Carlo simulations as £ 0.44 ± 0.01 
by Carlson et al. 19], and as £ = 0.42 ± 0.01 by As- 
trakharchik et al. |20j |. So far, simpler approaches have 
failed to agree on the value of £, and they deviate con- 



siderably from the Monte Carlo results. Heiselberg j2J 
obtained £ = 0.326, while Baker [H found £ = 0.326 and 
£ = 0.568 from different Pade approximations to Fermi 
gas expansions. Engelbrecht et al. |23| obtained £ = 0.59 
in a calculation based on BCS theory. 

It is the purpose of this paper to present a simple cal- 
culation that determines the universal constant £. It is 
based on Kohn-Sham density-functional theory (DFT) 
with a two parameter density functional that is fit to one 
analytically known result. The resulting value £ « 0.42 
is close to recent Monte Carlo results. This paper is or- 
ganized as follows. First, we present the DFT for the 
Fermi systems in the unitary regime. The density func- 
tional has a particularly simple and constrained form in 
the unitary regime. Second, we test and validate the 
density functional through calculations for the exactly 
solvable Calogero model. 

Carlson et al. El performed quantum Monte Carlo 
simulations for systems of N fermions in the unitary 
regime, with the number of fermions N ranging from 
N = 10 to approximately N s» 40. They found in par- 
ticular that the relation holds with very good accu- 
racy for all even number systems. This is a remarkable 
finding, since exact quantum mechanical energies usually 
differ from the corresponding Thomas-Fermi energies due 
to finite size effects and shell effects, both of which are 
apparently very small for fermions in the unitary regime. 
This suggests that the density functional from Thomas- 
Fermi theory 



(2) 



with c = jq (37r 2 ) 2 / 3 is a good approximation of the exact 
density functional, and that corrections might easily be 
accounted for via full-fledged Kohn-Sham DFT. 

As a starting point, we thus consider Thomas-Fermi 
theory of harmonically trapped fermions in the unitary 
regime. The density functional is 



cHO 
C TF 



[p] =£ TF [p} + ^mcj 2 r 2 p. 



(3) 



This functional is minimized under the condition that 
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the density is normalized to ./V particles. This yields the 
Thomas-Fermi density 



Ptf{t) 
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where Z = {Ti/muj) 1 / 2 is the oscillator length of the har- 
monic trap. For the Thomas-Fermi energy we insert the 
density (@J into the functional J2J and integrate over 
space. This yields 



E 



TF 



(5) 



For the two-fermion system, the exact quantum me- 
chanical result for the energy is 



E, 



2hu. 



(6) 



Let us assume that the relation also holds for harmon- 
ically trapped systems. Thus, we might equate Eq. © 
with Eq. © for N — 2, and solve for the universal con- 
stant £. This yields £ = 64/6 8/3 w 0.54. Note that this 
simple result deviates only about 20% from the Monte 
Carlo results 0, |2(j . This is quite encouraging and mo- 
tivates us to compute a more accurate estimate for £ via 
Kohn-Sham DFT. 

In Kohn-Sham DFT [25J, the ground-state density and 
energy of an interacting iV-fermion system is obtained 
from varying the (generally nonlocal and unknown) den- 
sity functional. Unfortunately, there is no simple recipe 
that permits one to construct the density functional. For 
dilute systems with sufficiently small range and positive 
scattering length, a systematic approach has been given 
by Puglia et al. [2(|. For electronic systems, one usu- 
ally parametrizes the density functional in terms of local 
densities and their gradients, and fits the considerable 
number of parameters to experimental data and theoret- 
ical results for infinite systems. This elaborate and cum- 
bersome approach has been successfully implemented in 
recent decades in quantum chemistry (see, e.g., Ref. [27J 
and references therein), and a similar approach is also 
pursued in nuclear structure |23, |2£J yfl Hl| . Note that 
DFT has also been successfully applied to study the BCS- 
BEC cross-over. In their study, Kim and Zubarev Q 
employed a density functional that was fit to known re- 
sults in the regime of very small scattering length and 
to the Monte Carlo results for the unitary regime, and 
employed a Pade approximation for intermediate values 
of the scattering length. 

Fortunately, the case of the unitary regime is much 
simpler. In what follows we restrict ourselves to small 
even-number systems with equal number of fermions in 
the two spin states. As the interaction does not in- 
troduce any new length scale into the system, a lo- 
cal density functional can only contain density terms 
proportional to p 5 ^ 3 h 2 /m, and gradients of the form 
p-(2k-2)/3^Q2kp-^ft w ^ n i n t e g er k > 0. Here, we al- 
low only for the density dependent term and incorporate 



gradient terms through an effective mass. Our ansatz for 
the density functional thus becomes 



( 4 ) £[P] 



m 



m 



2m 



eff 



N 

E 



cp 



5/3 



1^1 1 5 



(7) 

The univcr- 



The density is given as p = J2j 

sal constant £ and the effective mass m e g are parameters 
that will be determined below. The effective mass is, 
in principle, iV-dependent, but we suppress this depen- 
dency here. Note that this density functional has two 
important properties: First, its Thomas-Fermi limit is 
given in Eq. J5J) and is thus proportional to the density 
functional of the free Fermi gas. Second, non-localities 
of the density are introduced through the effective mass. 
Note that more elaborate approaches for superfluid sys- 
tems also introduce pairing densities [3^, IU 03 HH in 
the density functional in order to approximate the (un- 
known) non-local functional for superfluid systems. The 
quality of the results presented below, however, suggests 
that this is not necessary for the purpose of this study. 

For a determination of the parameters £ and m e ff we 
consider the problem of two spin-1/2 fermions inside a 
harmonic trap that interact via a zero-range interaction 
with infinite scattering length. The ground state of the 
two-fermion system is a spin singlet, and the relative co- 
ordinate wave function has been given by Busch et al. 

M 



A 



^cx(r) = (27r\ 3 y 1 ^-e- < - r /^ ' 2 . 



(8) 



Here, A = \ 2l is given in terms of the oscillator length I. 
The wave function diverges like 1/r for small distances 
due to the infinite scattering length; in practice, this di- 
vergence could be cut off by any non-zero range of the 
interaction potential (see, e.g. the discussion in Ref. 19] ), 
and it does not cause any problems in the normalization 
of the wave function. We recall that the ground-state en- 
ergy © of the two-fermion system is considerably lower 
than for nonintcracting fermions. Employing the relative 
wave function (JSJ) and the Gaussian ground state for the 
center of mass coordinate, we arrive at the density 



r 3/2 Z 3 



l _ e -2{r/lf 



dx t 



(9) 



A plot of this density is shown in Fig. E]as a full line, 
and can be compared to the noninteracting case (dotted 
line) . 

We add the term (m / 2)u> 2 r 2 p(r) of the harmonic con- 
finement to the density functional 10) and use Kohn- 
Sham theory to compute the density for given sets of 
parameters £ and m e ff. The best agreement with the 
exact density 10 and exact energy (JSJ is obtained for 
£ = 0.42 and m/m c ff = 0.69. The resulting density is 
plotted as the long-dashed line in Fig. ^ and exhibits 
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0.8 



exact 



^=0.44, m/m eff =0.52 
^=0.40, m/m, f =1 
£=0.42, mJm =0.69 

~ eff 




FIG. 1: (color online) Density of the harmonically trapped 
two-fermion system. Exact result (full line) compared to 
Kohn-Sham results from three different density functionals 
(dashed lines and dashed-dotted line), and the noninteract- 
ing system (dotted line). 



© yields £ = 4/9. Note that the resulting density is 
very close to the short-dashed line in Fig. ^ (labeled as 
i = 0.44, m/m cS = 0.52). 

In the unitary regime, the rather simple density func- 
tional J7J) does yield much improved and reliable results 
compared to Thomas-Fermi theory. To further test the 
form of this density functional we consider another inter- 
acting iV-body system whose density functional is also 
proportional to that of a free Fermi gas. The Calogero 
model 13611 with Hamiltonian 



,2 N 



H 



E 



1 N 

1 2 \ " 2 

-raw 2^x jt 



E 



§(§-!)' 



i d 2 



2 dx 2 f-f (xi - xj) 2 



(10) 



is exactly solvable in one dimension. Here, j3 > 1 is a 
dimensionless coupling constant. The exact ground-state 
energy of this Hamiltonian is 



only small deviations from the exact density. The energy 
deviates about 0.1% from its exact value @). Note that 
our value of the universal constant is very close to the 
Monte Carlo results [Til^. 

We want to estimate the robustness and stability of 
our result. For this purpose we consider the following 
two cases. First, we fixed the value of the universal con- 
stant to the Monte Carlo result £ = 0.44 by Carlson et al. 
|l9j . and obtain the effective mass m/m c s = 0.52 from 
the best fit to the exact density and energy. The resulting 
density is shown as the short-dashed line in Fig. ^ an d 
the energy deviates less than 1% from the exact result. 
Second, we fix the effective mass to m/m e s = 1, and ob- 
tain £ = 0.40 from the fit of the density functional. The 
resulting density is depicted as the dashed-dotted line in 
Fig.^ while the energy again deviates less than 1% from 
the exact result. These results show that DFT yields 
quite robust results for the universal constant £, and sys- 
tematically improves upon the naive Thomas-Fermi re- 
sult. This suggests that the good agreement of our best 
value £ = 0.42 with the Monte Carlo results is not merely 
accidental but due to the quality of the density functional 
we employed. 

We might even take our approach one step further and 
obtain a simple analytical estimate for the universal con- 
stant £. Due to the relatively large value of the effec- 
tive mass, the prefactor of the density-dependent term in 
the density functional (J7J becomes rather small, and we 
might neglect it by simply setting m/m c ff = £. In pres- 
ence of the confining harmonic potential, the Kohn-Sham 
equation is then identical to the Schrodinger equation of a 
three-dimensional harmonic oscillator where the kinetic 
energy is modified by a factor £. The analytical result 
for the ground-state energy of the two-fermion system is 
then E = STiu^ 1 / 2 , and comparison with the exact result 



' N 8 

E ex = hw[- + ^n(n - r 



(11) 



Let us focus on the two-body interaction. The inverse 
square potential does not introduce any new length scale 
into the Hamiltonian as it scales as the kinetic energy. 
Thus, the inter-particle distance is the only length scale, 
and the density-dependent energy must be proportional 
to the one-dimensional Fermi gas, the proportionality 
constant being denoted as i] 2 



£tf[p] 
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(12) 



In this respect, the Calogero model is similar to the Fermi 
gas in the unitary regime. This approach leads dir ectly 
to the Thomas- Fermi theory for the Calogero model [37| • 
The constant rj 2 can be determined from Thomas- 
Fermi theory once we add the confining harmonic po- 
tential and thereby make contact with exactly known re- 
sults. Thus, the density functional becomes 



£[p] = £ T f[p\ + 7; muj2 x 2 p, 
and it is minimized by the normalized density 



Ptf 



(r/lf 



(13) 



(14) 



Here, I = {Ti/muj) 1 / 2 again denotes the oscillator length. 
The Thomas-Fermi density l|14f) is Wigner's semi circle, 
and agrees with the standard approach pjsj . 

We insert the density l|14f> into the functional (|13|l and 
integrate. The resulting Thomas-Fermi energy is 



-HluN 2 



(15) 



4 



and comparison with the exact result Ijllfl fixes T) = (3/2. 
Note that the Thomas-Fermi energy (II 51) differs from the 
exact result l|ll|) by considerable finite size corrections, 
and this is a difference to the Fermi gas in the unitary 
regime. 

Let us apply Kohn-Sham DFT to the Calogero model. 
The ansatz for the density functional is in full analogy to 
the one we made for the Fermi gas in the unitary regime 



m = h - 

m 



o^ei^)i 2 +^G 2 -— y 

2m off f-f 6 V m cff / 



+ —muj 2 x 2 p(x) . 



(16) 



Here we have again introduced the effective mass m e s as 
the only fit parameter and have already added the term 
due to the harmonic confinement. 

We determine the effective mass by solving the Kohn- 
Sham equation for the density functional i|16|) and com- 
pare the resulting density with exact results. Note that 
the exact density of the Calogero model (|10|l is only 
known for a few values of the coupling constant /3, though 
the many-body ground state of this model is known for 
decades. For (3 = 1, /3 = 2, and (3 — 4, the density 
is related to the eigenvalue distribution of the orthogo- 
nal, unitary and symplectic Gaussian random matrix en- 
semble, respectively. Analytical expressions are given in 
Refs.|39ll40l|. We focus on the case (3 = 4, and determine 
the effective mass by fit. From calculations for particle 
numbers N — 2, 4, 16, 32 one obtains approximately 



6.3 + 



8J3 



(17) 



Thus, the effective mass is considerably smaller than the 
mass. Figure El shows that the Kohn-Sham densities are 
close to the exact results. The Thomas-Fermi result is 
also shown for comparison. Note that only a relatively 
small effective mass (II Til reproduces the density oscilla- 
tions. Note also that the deviation of the DFT energies 
from the exact result is about half as large as the error 
of the corresponding Thomas-Fermi energies. This shows 
that the simple density functional (|16|) yields significantly 
improved energies and densities. 



In summary, we have used density-functional theory 
to compute the universal constant of the Fermi gas in 
the unitary regime. This approach is based on the ob- 
servation that the Thomas-Fermi energy is a reasonable 
first order approximation to the quantum mechanical re- 
sults, and on the constraints that the unitary regime 
imposes on the form of the density functional. Our es- 
timate £ = 0.42 results from a best fit to the density 
and energy of the harmonically trapped two-fermion sys- 
tem, and is in good agreement with much more elaborate 
Monte Carlo studies. Our result is stable with respect to 
variations of the density functional, and favors a suffi- 
ciently large effective mass. The particular form of the 
density functional could also be tested in applications to 
the Calogero model. 
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FIG. 2: (color online) Densities of the Calogero model for 
(3 — 4 for different numbers of particles N. Full line: density 
functional theory; long-dashed line: leading order expansion 
of exact density; short-dashed line: Thomas-Fermi theory. 
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